Can the SO (10) Model with Two Higgs Doublets 



Reproduce the Observed Fermion Masses? 
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It is usually considered that the SO(10) model with one 10 and one 126 Higgs scalars cannot 
reproduce the observed quark and charged lepton masses. Against this conventional conjecture, we 
find solutions of the parameters which can give the observed fermion mass spectra. The SO(10) 
model with one 10 and one 120 Higgs scalars is also discussed. 
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I. INTRODUCTION 

The grand unification theory (GUT) is very attractive 
as a unified description of the fundamental forces in the 
nature. Especially, the SO(10) model is the most attrac- 
tive to us when we take the unification of the quarks and 
leptons into consideration. However, in order to repro- 
duce the observed quark and lepton masses and mixings, 
usually, a lot of Higgs scalars are brought into the model. 
We think that the nature is simple. What is of the great- 
est interest to us is to know the minimum number of the 
Higgs scalars which can give the observed fermion mass 
spectra. A model with one Higgs scalar is obviously ruled 
out for the description of the realistic quark and lepton 
mass spectra. Then, how is a model with two different 
types of Higgs scalars (e.g., 10 and 126 scalars)? 

In the SO(10) GUT scenario, a model with one 10 and 
one 126 Higgs scalars leads to the relation |jj 



M e = c u M u + c d M dl 



(1.1) 



where M e , M u and M d are charged lepton, up-quark and 
down-quark mass matrices, respectively. It is widely ac- 
cepted that there will be almost no solution of c u and 
c d which give the observed fermion mass spectra. The 
reason is as follows: We take a basis on which the up- 
quark mass mat rix M u is diagonal (M u — D u ). Then, 
the relation ([hi]) is expressed as 



M e = c u D u + c d M d . 



(1.2) 



Considering that M d is almost diagonal and the mass hi- 
erarchy of up-quark sector is much severe than that of 
down-quark sector, we observe that the contribution to 
the first and the second generation part of M e from the 
up-quark part D u is negligible so that it is proportional 
to that of M d . Thus, the relation (1.1) which predicts 
m e /m^m d /m s does not reproduce the observed hier- 
archical structure of the down-quark and charged lepton 
masses such as predicted by Georgi- Jarlskog mass re- 
lations rrift = m T , m s = m M /3 and m d = 3m e at the 



GUT scale ||. However, the above conclusion is some- 
what impatient one. (i) It is too simplified to regard M d 
as almost diagonal, (ii) We must check a possibility that 
the mass relations are satisfied with the opposite signs, 
i.e., = ±m r , m s = and m d — ±3m e . (iii) 

The mass values at the GUT scale which are evaluated 
from the observed values by using the renormalization 
group equations show sizable deviations from the Georgi- 
Jarskog relations. The purpose of the present paper is to 
investigate systematically whether there are solutions of 
c u and c d which give the realistic quark and lepton masses 
or not. 



II. OUTLINE OF THE INVESTIGATION 

In the SO(10) GUT model with one 10 and one 126 
Higgs scalars, the down-quark and down-lepton mass ma- 
trices M d and M e are given by 



Md = M + Mi, M e = M -3M u 



(2.1) 



where Mq and Mi are mass matrices which are generated 
by the 10 and 126 Higgs scalars 4>w an d </>i26, respec- 
tively. Inversely, we obtain 



M 



1 



1 



(3M d + M e ), Mi = -{M d -M e ). (2.2) 



On the other hand, the up-quark mass matrix M u is given 
by 



CqMo + CiMi, 



where 



co = </v d a = «)/(^ °), 

Ci = V^/vf = (<t>$ 6 ) 



(2.3) 



(2.4) 



and <p u and 4> d denote Higgs scalar components which 
couple with up- and down-quark se ctor s, respectively. 
Therefore, by using the relations Eq.(2.2), we obtain the 
relation 



1 



where 



Cd 



M e = c d M d + c u M u , 



3c + ci 



c - ci 



c - ci 



(2.5) 



(2.6) 



For convenience, first, we investigate the case that the 
matrices M u , M d and M e are symmetrical matrices at the 
unification scale because we assume that they are gener- 
ated by the 10 and 126 Higgs. Then, we can diagonalizc 
those by unitary matrices U u , U d and U e , respectively, as 



UlM u U u = D u , UjM d U d = D d , Ul M e U e 



(2.7) 



where D u , D d and D e are diagonal matrices. Since the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix V is given 
by 



v = K u* d , 



(2.8) 



the relation (pi]) is re- written as follows: 

{UtU u ) T D e {UtU u ) = c d VD d V T + c u D u . (2.9) 

At present, we have almost known the experimental val- 
ues of D e , D u and VD d V^. Therefore, we obtain the 
independent three equations: 



TrD e Dl = \c d \ 2 Tr[(VD d V T + nD u ){VD d V T + kD u )^ 
Ti(D e Dl) 2 = \c d \ i Ti[((VD d V T + KD a )(VD d V T + kD u )^ 2 
detD e Dl = \c d \ 6 det \{VD d V T + nD u ){VD d V T + kD u )^ 



(2.10) 
(2.11) 
(2.12) 



where k = c u /c d . By eliminating the parameter c d , we 
have two equations for the parameter k: 



{ml + ml + mlf (2.10) 



(m 2 . 



to' 



~ 2 2 2 

m^m^m^ 
+ to 2 ) 2 



J2.12| ) ' 
(f2~10l) 2 



2(to 2 to 



(I2.10D 2 - (|2.11|) 



(2.13) 



(2.14) 



where ( 2.1C ) 3 , for instance, means the right-hand side of 
(2.10) to the third power. Let us deno te th e parameter 
values of k evaluated from ( 2.13 ) and ( 2.14 ) as ka and 
Kb, respectively. If ka and kb coincide with each other, 
then we have a possibility that the SO(10) GUT model 
can reproduce the observed quark and lepton mass spec- 
tra. If ka and kb do not so, the SO(10) model with 
one 10 and one 126 Higgs scalars is ruled out, and we 
must bring more Higgs scalars into the model. Of course, 
in the numerical evaluation, the values ka and Kb will 
have sizable errors, because the observed values D e , D u , 
D d and V have experimental errors, and the values at 
the GUT scale also have errors. The values ka and kb 
are not so sensitive to the renormalization group equa- 
tion effect (evolution effect), because those are almost 
determined only by the mass ratios. (More details will 
be discussed in the Sec. III.) Therefore, we will evaluate 
ka and Kb by using the center values at /i = mz in the 
Sec. IV. If we find ~ kb, we will give further detailed 
numerical study only for the case. 



III. EVOLUTION EFFECT 

The relations (2.13) and (2.14) hold only at the unifica- 
tion scale /i = Ax On the other hand, we know only the 
experimental values of the fermion masses m f and CKM 
matrix parameters Vij at the electroweak scale [i = mz- 
For a model which does not have any intermediate en- 
ergy scales, we can straightforwardly estimate the values 
of m/ and Vij at = Ax from those at fi = mz by the 
one-loop renormalization equation 



dYj_ 

dt 



167T 



;(T f -G f + H f )Yf 



(3.1) 



where Tf,Gf and Hf denote contributions from fermion- 
loop corrections, vertex corrections due to the gauge 
bosons and vertex corrections due to the Higgs boson(s), 
respectively. Therefore, we can directly check the re- 
lations (2.13) and (2.14) by substituting the observable 
quantities mf and Vij at /i = Ax- However, for a model 
which has an intermediate energy scale such as a non- 
SUSY model, the values of m/ and Vij at fi — Ax are 
highly model-dependent, so that the check of Eqs. (2.13) 
and (2.14) cannot be done so straightforwardly. 

In this section, we will show that we can approximately 
check Eq. (2.13) and (2.14) by using the values of to/ and 
Vij at — mz, without knowing the explicit values of to/ 
and Vij at /i = Ax, as far as the evolutions of to./ and 
Vij are not singular. 

It is well known that in such a conventional model the 
evolution effects are approximately described as m 



mP u jrr^ t m c jm\ 



m u /m t m c /m t 



1 



2 



m s/ m b 



m d/ m b 
m d /m b 

\v ub \ \v cb \ \v td \ \v ts \ 

m°Jm° c „m° d /m° a ^ |V£| 



W c ° b \ 



l+e d 



\v t ° s \ 



k° = v°Dy oT + k u d: 



\ V o d \ „ 



m u /m c m d /m s \V US \ \V cd \ 



1 , 



where mP q and V-j (m q and Vij) denote the values at 
li = Ax (m = m z)- The relations (3.2) hold only for 
a model where the Yukawa coupling constant of top 
quark, y t = (Y„) 33 , satisfies y t > {Y d )ij = 1,2,3). 
The relations (3.2) also hold even in a model which 
has an intermediate energy scale A/, because, for exam- 
ple, when we denote (m u /mt) /J= A x /(m u /mt) /J =A J and 
(m u /m t ) M =A 7 /(m u /m t ) M=mz as l+e ul and l+e u2 , re- 
spectively, we can obtain {m u /m t )^=A x / 'K/mt),i=m z — 
1 + e u with e u = £ u i + e u2 . 

By using the approximate relations (3.2) the diagonal- 
ized up-quark mass matrix D® at /j, = Ax is presented 
as 



m^/m 



m+ 



ml 



mi 






m u /m t 







m® c jm\ 
1 



m c /m t 
1 




= -L(l + e u S)D U) 
rot 



where 




Similarly, the matrix D° is given by 

D° d ^^(l + s d S)D d . 
m b 

The CKM matrix V° at = A x is given by 



(3.3) 



(3.4) 



(3.5) 



1 V us V ub (l + e d ) 

V°~\ V cd 1 V cb (l+e d ) 

V td (l + e d ) V ts (l + e d ) 1 

~ (l + e d S 3 )V(l + e d S 3 ) -2e d S 3 , 



(3.6) 



where S 3 = 1 — S and 1 is a 3 x 3 unit matrix. By using 
the relations (3.4) - (3.6), we can obtain the approximate 
expression 



V°D° d V 0T 



m b 



[(1 + e d )VD d V T - e d m b S 3 ] , (3.7) 



where we have used the observed hierarchical relations 
among the quark mass ratios and CKM matrix parame- 
ters. Therefore, the matrix VD d V T + nD u in Eqs. (2.10)- 
(2.12) is given by 



OnO 



m 



o 



(l+e d )^(VD d V T + kD u 
m b y 

-e d m b S 3 + s u kD u S) , 



(3.2) where 



k = 



rrq/mt 



m\jm b 1 + e d 



(3.8) 



(3.9) 



Since the solutions n are of the order of 10~ 2 as we 
show in the next section, we can neglect the term kD u S 
compared with VD d V T (note that in order to neglect 
the component (D u S)n it is essential that the sign of 
m d /m s is positive, because (VD d V T )n ~ m d + V^ s m s 
and V^s — \ m d/m s \). On the other hand, for such a small 
value of k, the term m b S 3 cannot be neglected compared 
with the term kD u . However, for a small value of s d , 
we can find that the solutions k are substantially not af- 
fected by the term e d m b S 3 . As a result, we obtain the 
approximate expression 



(l + e d )^(VD d V T + kD u ). 
m b 



Therefore, Eq. (2.13) and (2.14) at (j, = A x , i.e., 
[K) 2 + (m°) 2 + K) 2 ] 3 [Tr(^°^°t])] 3 



«) 2 K) 2 K) 2 



[K) 2 + (m°) 2 + K) 



,0\2l2 



2[(m°) 2 (m°) 2 + (m») 2 K) 2 + K) 2 K) 2 ] 



(3.10) 



(3.11) 



(3.12) 



[Tr(K°K°t)] 2 - Tr(K°K°i) 2 ' 
are approximately replaced by the relations at [i = mz- 
(ml + ml + m 2 ) 3 _ [Tt(KK^)} 3 



mlm^ml 



det(KKi) 



(m 2 + 



^?) 2 



2(m 2 m 2 + m 2 m 2 + m 2 m 2 ) 



(3.13) 



(3.14) 



where 



Tr[(KKl)] 2 - Tr(KK^) 2 ' 



K = VD d V^ + kD v 



(3.15) 



and k is given by Eq. (3.9). This means that when we 
find the solution k at fj, = mz, the solution at \x = Ax 
also exists, no matter whether the model is a SUSY one 
or a non-SUSY one. Then, we can obtain the value k° 
at \jl = Ax from the relation (3.9) with the solution k at 
H = m Z - 
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IV. NUMERICAL STUDY AT p = Mz 

As mentioned in the preceding section, if the solution 
k exists at the energy scale p = mz, the one at p = Ax 
also exists. Therefore, we investigate th e rela tions (2.13 ) 
and (2.14) at p = m z . Note that Eqs.(fjf) and ( |2.14|) 
are realized by GUT scale because Eq.(2.7) is broken 
at p — mz- In the present section, tentatively, we as- 
sume that the Yukawa coupling constant Yio and 3^26 a t 
p = mz keep their forms symmetrical, so that we can 
put the obser ved v alues D u , D g and V at p, = m z into 
the relations ( [2.13 ) and (2.14). For the fermion masses 
at p — mz, we use the following values: H 



677t^ MeV, 
2.33±°;g MeV, 
1746.7 ±0.3 MeV, 



m t = 181 ± 13 GeV, 



m b = 3.00 ±0.11 GeV, 



m d 



93.4±&g MeV, 
4.69™ MeV, (4.1) 



to m = 102.75138 ±0.00033 MeV, 
m e = 0.48684727 ± 0.00000014 MeV. 



The input values for the CKM matrix parameters have 
been taken as [p| 



0l2 
9l3 



0.219 
0.002 



0.226, 
0.005, 



^23 



0.037- 0.043, 



(4.2) 



where 



C13C12 C13S12 Si3e 



V = -C23S12 - S23Cl2Sl3e id C23C12 - S23Sl2Sl3e ^ ' :, S23C13 , (4.3) 



S23S12 - C23Cl2Sl3e* 5 -S23C12 - C23Sl2Sl2Sl3e l5 C23C13 



with Cij = cos 9ij and Sjj = sin Oij . The calculation 
has been performed allowing all the combinations of the 
quark mass signatures. Here it should be noted that, 
since m u is much smaller than m c and mt , the difference 
of the sign of m u scarcely makes a change of allowed re- 
gions. In this calculation, we have selected $23 and 5 as 
input parameters and m s , Cd and k as output parameters 
because the calculation is sensitive to these parameters. 
We give the numerical results in Fig 1. Here, except 
for to.,, #23 and 5, we have adopted the center values of 



Eq.(4.1) as input values. Moving $23 at intervals of 0.0005 
rad and fixing S — 60°, we search the solutions where ka 
and Kb become coincident. Our numerical analysis shows 
that the solutions exist in the combinations of Table Q. 
In a table [n| we show the nearest solution of m s , 623 and 



S to the center values of Eq.(4.1). 

In the following we perfor m dat a fittin g for the case 
of top line of Table || Eqs. j2~To|)-(|2~ll) can constrain 
only the absolute value of c^. The argument of the pa- 
rameter c ( i may be decided by taking neutrino sector into 
consideration in the future. For the time being, we set 
so that cq becomes a real number: 



Cd 



\Cd\e 



ia _ „0.107i 



1 — Cd 

c = = 34.7, 



ci = 



Cu 

3 + c d 



= 101.8 — 10.8*. 



(4.4) 
(4.5) 



In this case, the mass matrices in MeV are 



01/ri t/T 1 , n 1 T/n \rT\ / -12.4 - 0.7« -23.0 - \M 9.6 - 13.2* 

Mo = WDdV + Cd{ f u + VDdV } = I -23.0 - 1.8i -91.5 - 3.9* 194.0 + 10.5* ) . (4.6) 

9.6-13.2*' 194.0 + 10.5* 1874.9- 180.0* 

vn ,-, f) , /) ! ' / -*-19 + 0.69i 7.68 + 1.43* -3.72 + 4.09* 

M x = -=^ c d {KiJ u +vUdV ) = j 7M + x A3i 24.14 + 3.88* -65.05 - 10.48* |. (4.7) 

4.09*' -65.05- 10.48* 1119.67+ 179.98*' 




Here, using the condition \Z\vq \ 2 -i 
246GeV, we can get VEV's as 

246 [GeV] 



-,d\2 



„«|2 



V(|c | 2 + l) + (|ci| 2 + l)|p| 5 



(4.8) 



with p = vf/vg. Then, the Yukawa couplings about 10 
and 126 become 



M 



Mi 



126 



(4.9) 



We consider that the model should be calculable pertur- 
bativly. We can see that every element of the Yukawa 
coupling constants (4.9) is smaller than one if we take a 
suitable value of \p\. 
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V. 10 AND 120 



VI. SUMMARY AND DISCUSSION 



In the SO(10) GUT scenario, we can also discuss the 
model with one 10 and one 120 by the same method. 
The Yukawa couplings of 10 and 120 are symmetric 
and antisymmetric, respectively. If we consider a case 
that the Yukawa coupling constants of 10 are real and 
120 pure imaginary, we can make them Hermitian, i.e., 
Yjo = Yio and Y^ 2a — Yi 2 q. Therefore, by considering 
the real vacuum expectation values vio and V120, we can 
obtain the Hermitian mass matrices M u , Md and M e : 

M d = Mo + M 2 , M e = M - 3M 2 , 

M u = c M + c 2 M 2 . (5.1) 

Then, we can diagonalize those by unitary matrices U u , 
Ud and U e as 

UtM u U u = D u , U\M d U d = D d , U\M e U e = D e . 

(5.2) 

Since the CKM matrix V is given by 

V = UtU d , (5.3) 

the relation (|]l]) is re- written as follows: 

(UlU e )D e (UtU e )^ = CdVDdV^ + c u D u . (5.4) 

As stated previously, we have almost known the exper- 
imental values of D e , D u and VDdV^ . Therefore, we 
obtain the independent three equations: 



TrD e = c d [TiD d + kTtD u ] 



(5.5) 



TiD 2 e = c d [TrDj + 2 K Ti{D u VD d V^) + k 2 TiD u ], (5.6) 



det£> e = c 3 d det(VD d V^ + kD u ) 



(5.7) 



where k = c u /c d - For the parameter k, we have two 
equations: 



nil + m fi 



(m e + m M + m T ) 2 



_ Tr£> 2 + 2/sTr(A t U£> d I/t) + K 2 Tr£> 2 
~ (TrD d + KTvD u y : 

m e m u m T det(VD d V^ + kD u ) 



m. 



m T ) 3 (TtDd + KTrD u y 



(5.8) 
(5.9) 



Eqs. (|5.8| ) and ( |5.9|) are more simple than Eqs. ( 2.13| ) 
and (2.14). Cd and n are real since we have assumed the 
M u , Md and M e to be Hermitian. So the calculation 
is easier than the case for 10 and 126. The numerical 
results are listed in Table nil 



In conclusion, we have investigated whether an SO(10) 
model with two Higgs scalars can reproduce the observed 
mass spectra of the up- and down-quark sectors and 
charged lepton sector or not. What is of great inter- 
est is to see whether we can find reasonable values of 
the parameters c u and Cd which satisfy the SO(10) rela- 
tion (2.5) or not. For the case with one 10 and one 126 
scalars, in a par amete r k = c u /cd, we have obtained two 
equations (2.13) and ( p. 14 ) which hold at the unification 
scale fi = Ax and which arc described in terms of the 
observable quantities (the fermion masses and CKM ma- 
trix parameters). We have sought for the solution of k 
approximately by using the observed fermion masses and 
CKM matrix parameters at fi = mz instead of the ob- 
servable quantities at u, = Ax- Although we have found 
no solution for real k, we have foun d four solut ions for 
complex k which satisfy Eqs. ( 2.13| ) and ( p. 14 ) within 
the experimental errors. Similarly, we have found four 
solutions for a model with one 10 and one 120 scalars. 
It should be worth while noting that the solutions in the 
latter model are real. The latter model is very attractive 
because the origin of the CP violation attributes only to 
the 120 scalar. In the both models, we can make the 
magnitudes of all the Yukawa coupling constants smaller 
than one, so that the models are safely calculable under 
the perturbation theory. 

By the way, note that the numerical results are very 
sensitive to the values of m s and #23- For numerical fit- 
tings, it is favor that the strange quark mass m s is some- 
what smaller than the center value m s = 93.4 MeV which 
is quoted in Ref. ||. 

Also note that the relative sign of md to m s in each 
solu tion is positive, i.e, md/m s > as seen in Tables || 
and III. It is well known that a model with a texture 



(Md)n = on the nearly diagonal basis of the up-quark 
mass matrix M u leads to the relation | V us \ = *J —md/m s 
[Q, where the relative sign is negative, i.e., mdjm s < 0. 
On the contrary, we can conclude that in the SO(10) 
model with two Higgs scalars, we cannot adopt a model 
with the texture (Md)n = 0. 

In the present paper, we have demonstrated that 
the unified description of the quark and charged lep- 
ton masses in the SO (10) model with two Higgs scalars 
is possible. However, we have not referred to the neu- 
trino masses. Concerning this problem, Brahmachari 
and Mohapatra have recently showed that one 10 and 
one 126 model is incompatible with large v u -v T mixing 
angle Q. Since there are many possibilities for neutrino 
mass generation mechanism, we are optimistic about this 
problem, too. Investigating for a question whether an 
SO(10) model with two Higgs scalars can give a unified 
description of quark and lepton masses including neu- 
trino masses and mixings or not is our next big task. 
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num. 


(mt 


m c ,m u ) 


(m b ,m a ,m d ) 


(m T ,m M ,m e ) 


(a) 


(+ 


- +) 


(+ - -) 


(+ ± ±) 


(b) 


(+ 


- -) 


(+ - -) 


(+ ± ±) 



TABLE I. The combinations of the signs of (m t ,m c ,m u ), 
(mb,m a ,md) and (m T ,m M ,m e ). The notation (mt, m C ] m u ) 
= (+-_+) denotes m t > 0, m c < and m u > 0. Eqs. (2.12) 



and (2.14) are not affected by the signs of charged leptons. 





Input 






Output 






|ft»|[rad] 


S[°] 


m s [MeV] 


Cd K 




(a) 


0.0420 


60. 


76.3 


3.15698 -0.01928 - 


0.00089i 




0.0420 


60. 


76.3 


3.03577 -0.01937 - 


O.OOlOli 


(b) 


0.0420 


60. 


76.3 


3.13307 -0.01929 - 


0.00092i 




0.0420 


60. 


76.3 


3.00558 -0.01939 - 


0.00105i 



TABLE II. Four sets of parameters giving good data fit- 
ting at fi = mz for one 10 and one 126 Higgs scalars. (a) 
and (b) correspond to the mass signatures in Table I, and the 
upper and lower lines do to the two intersections in Fig.l 



num. 


(mt 


m c ,m u ) 


(m b ,m s ,m d ) 


(m T ,m M ,m e ) 


(a-1) 


(+ 


- +) 


(+ - -) 


(+ + +) 


(a-2) 


(+ 


- +) 


(+ - -) 


(+ + -) 


(b-1) 


(+ 


" ") 


(+ - -) 


(+ + +) 


(b-2) 


(+ 


- -) 


(+ - -) 


(+ + -) 



TABLE III. The combinations of the signs of (ra t , m c , m u ), 
(nib, m s ,md) and (m T , m M , m e ) for one 10 and one 120 Higgs 
scalars. 





Input 






Output 






|0 23 |[rad] 


S[°] 


m s [MeV] 


Cd 


K 


(a-1) 


0.0415 


60. 


79.551 


0.05905 


— 0.01957 


(a-2) 


0.0415 


60. 


79.238 


0.06124 


— 0.01942 


(b-1) 


0.0415 


60. 


79.673 


0.05855 


-0.01960 


(b-2) 


0.0415 


60. 


79.316 


0.06080 


-0.01945 



TABLE IV. Four sets of parameters giving good data fit- 
ting at n = mz for one 10 and one 120 Higgs scalars. (a-i) 
and (b-i) correspond to the mass signatures in Table III. 



0.0015 



0.001 



0.0005 - 



- 



-0.0005 



-0.001 



-0.0015 



I * 1 1 1 I 1 1 1 1 I 1 1 1 1 I 



-0.020 



-0.019 -0.018 

Re K 



-0.017 



FIG. 1. The relations between Eqs. (2. 13) and (2.14) on the 
complex plane of ft. The solid (dotted) line shows the solution 
of Eq.(2.13) (Eq.(2.14)). 
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